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5— i " Abstract 
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We covariantize the decoupling limit of massive gravity proposed in [1] and 
study the cosmology of this theory as a proxy, which embodies key features 
of the fully non-linear covariant theory. We first confirm that it exhibits a 
self-accelerating solution, similar to what has been found in [2], where the 
Hubble parameter corresponds to the graviton mass. For a certain range of 
Q_[ parameters fluctuations relative to the self-accelerating background are stable 

and form an attractor solution. We also show that a degravitating solution 
can not be constructed in this covariantized proxy theory in a meaningful 
way. As for cosmic structure formation, we find that the helicity-0 mode of 
■ the graviton causes an enhancement relative to ACDM. For consistency we 

also compare proxy theories obtained starting from different frames in the 
^ [ decoupling limit and discuss the possibility of obtaining a non-representative 

£T) ' proxy theory by choosing the wrong starting frame. 
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jxj ■ 1 Introduction 

Observations of the CMB, supernovae, lensing and baryon acoustic oscillations have 
led to the cosmological standard model which requires an accelerated expansion of 
the late Universe, driven by dark energy but despite many years of research its origin 
has not yet been identified. There are two major explanations for the origin and 
properties of dark energy. 

The first solution consists of introducing a cosmological constant A with a con- 
stant energy density causing an effective repulsive force between cosmological objects 
at large distances. From the particle physics perspective the cosmological constant 
could correspond to the vacuum energy density. The theoretical expectations for 
the vacuum energy density caused by fluctuating quantum fields, however, exceeds 
the observational bounds on A by up to 120 orders of magnitude. This discrepancy 
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remains for more almost a century one of the most challenging puzzles in physics, 
[3]. 

Alternatively, the acceleration of the Universe can be explained by introducing 
new dynamical degrees of freedom, either by invoking new fluids with negative 
pressure or by changing the geometrical part of Einstein's equations. In particular, 
weakening gravity on cosmological scales could not only tackle the cosmological 
constant problem, but would also come hand in hand with new degrees of freedom 
which might be responsible for a late-time speed-up of the Hubble expansion. Such 
scenarios could arise in massive gravity or in higher-dimensional frameworks. 

In the higher dimensional picture, the Dvali-Gabadadze-Porrati (DGP) model is 
one of the important large scale modified theories of gravity [4]. In this braneworld 
model our Universe is confined to a three-brane embedded in a five-dimensional bulk. 
On small scales, four- dimensional gravity is recovered due to an intrinsic Einstein 
Hilbert term sourced by the brane curvature, whereas on larger, cosmological scales 
gravity is systematically weaker as the graviton acquires a soft mass m which limits 
its effective range. Being a fundamentally higher dimensional theory, the effective 
four-dimensional graviton on the brane carries five degrees of freedom, namely the 
usual helicity-2 modes, two helicity-1 modes and one helicity-0 mode. Whilst the 
helicity-1 mods typically decouple, the helicity-0 one can mediate an extra fifth 
force. In the limit m — > 0, one recovers General Relativity (GR) through the 
Vainshtein mechanism: The basic idea is to decouple the additional modes from 
the gravitational dynamics via nonlinear interactions of the helicity-0 mode of the 
graviton, [5]. As a result, at the vicinity of matter, the non-linear interaction for 
the helicity-0 mode become large and hence suppresses its coupling to matter. This 
decoupling of the nonlinear helicity-0 mode is manifest in the limit where M4, M5 — > 
00 and m — > while the strong coupling scale A = (Mpim 2 ) 1 / 3 is kept fixed. This 
limit enables a linear treatment of the usual helicity-2 mode of gravity while the 
helicity-0 mode it is described non-linearly, which is the so-called decoupling limit. 

One of the successes of the DGP model is the existence of a self-accelerating 
solution, where the acceleration of the Universe is sourced by the graviton own 
degrees of freedom (more precisely its helicity-0 mode). Unfortunately that branch 
of solution seems to be plagued by ghost-like instabilities [6, 7, 8], in the DGP model, 
but this issue could be avoided in more sophisticated setups, for instance including 
Gauss-Bonnet terms in the bulk [9]. 

More recently, it has been shown that the decoupling limit of DGP could be 
extended to more general Galilean invariant interactions [10]. This Galileon model 
relies strongly on the symmetry of the helicity-0 mode tt: Invariance under internal 
Galilean and shift transformations, which in induced gravity braneworld models can 
be regarded as residuals of the 5-dimensional Poincare invariance. These symmetries 
and the postulate of ghost-absence restrict the construction of the effective 7r La- 
grangian. There exist only five derivative interactions which fulfill these conditions. 
From the five dimensional point of view these Galilean invariant interactions are 
consequences of Lovelock invariants in the bulk of generalized braneworld models, 
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[11]. Since their inception there has been a flurry of investigations related to self- 
accelerating de Sitter solutions without ghosts [10, 12], Galileon cosmology and its 
observations [13, 14], inflation [15, 16, 17, 18], lensing [19], superluminalities arising 
in spherically symmetric solutions around compact sources [20], K-mouflage [21], 
Kinetic Gravity Braiding [22], etc... . Furthermore, there has been some effort in 
generalizing the Galileon to a non-flat background. The first attempt was then to 
covariantize directly the decoupling limit and to study its resulting cosmology [13]. 
In particular it was shown in [23] that the naive covariantization would yield ghost- 
like terms at the level of equation of motion but a given unique nonminimal coupling 
between tt and the curvature can remove these terms resulting in second order of 
equations of motion [23], which are also consistent with a higher- dimensional con- 
struction [11]. In this paper, we will pursue the same strategy when constructing 
our proxy theory. The outcome of this kind of covariantization method was explored 
further in [24]. While this covariantization is ghost-free, the Galileon symmetry is 
broken explicitly in curved backgrounds. Only recently there has been a success- 
ful generalization to the (Anti-) de Sitter background and ultimately to maximally 
symmetric backgrounds where it has been discovered that the de Sitter Galileon in- 
teractions acquire additional, potential-like terms being functions of tt and mixtures 
of 7r and gradients of tt, all fulfilling a generalized Galileon symmetry [25]. 

There exists a parallel to theories centered on a massive graviton: Galileon-type 
interaction terms naturally arise in gravitational theories using a massive spin-2 
particle as an exchange particle, which has, in addition, been constructed to be 
ghost-free be it in three dimensions, [26] or for an generalized Fierz-Pauli action in 
four dimensions [1, 27] 3 . Not only is the existence of a graviton mass a fundamental 
question from a theoretical perspective, it could also have important consequences 
both in cosmology and in solar system physics, [32, 33]. Although solar system ob- 
servations have confirmed GR to high accuracy and placed bounds on the graviton 
mass to be smaller than a few ~ 10~ 32 eV, even such a small mass would become rel- 
evant at the Hubble scale which corresponds to the graviton Compton wavelength. 
In particular, it has been successfully shown that this massive gravity theory ex- 
hibits a stable self-accelerated solution in the decoupling limit since the scalar mode 
can generate a constant negative pressure density. In the decoupling limit, the 
expansion history of the Universe in this self-accelerating branch was found to be 
indistinguishable from ACDM, [2]. 

While the self-accelerating solutions in the above models yield viable expansion 
histories including late-time acceleration, they do not address the cosmological con- 
stant problem, i.e. the giant mismatch between the theoretically computed high 
energy density of the vacuum and the low observed value. A possible answer comes 
from the idea of degravitation, which asserts that the energy density could be as 
large as the theoretically expected value, but would not bear a large effect on the 

3 Such a theory was also constructed using auxiliary extra dimensions, [28, 29]. While in its 
most fundamental form, a ghost appears at quartic order in the decoupling limit, [30], it can also 
be cured order by order, [31]. 
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geometry. Technically, gravity is less strong on large scales (IR-limit) and could act 
as a high-pass filter suppressing the gravitational effect of a potentially large vac- 
uum energy. Since such modifications of gravity in the IR naturally arise in models 
of massive gravity, they logically provide a possible mechanism to degravitate the 
vacuum energy density, [34, 35, 36], which was observed in bi-Galileon models [37] 
as well as in the decoupling limit of massive gravity, [2]. Analogously, the DGP 
braneworld model can be extended to higher dimensions to tackle the cosmological 
constant problem as well, [36, 38, 39]. 

In this paper, we focus on the covariantization of a ghost-less extension of Fierz- 
Pauli massive gravity recently proposed in [27] and show that this proxy model 
allows for a stable self-accelerating solution. Hereby, we have performed the covari- 
antization in the Jordan as well as in the Einstein frame. We discuss the differences 
between the two approaches and the consequences of choosing the wrong starting 
frame. In the well defined proxy theory we are able to tackle the puzzle about 
the self-acceleration of the Universe but not the one about the cosmological con- 
stant problem. Furthermore, we study the perturbations around the self-accelerating 
background and provide the expression for the modified evolution equation for these 
density perturbations. As expected, the ir field enhances the gravitational clustering 
resulting in a rapid growth of structures. This result is quite different from the one 
we had obtained in [2], where the self-accelerating solution was indistinguishable 
from a ACDM, but we emphasize that the theory used in this paper, is a proxy 
model distinct from massive gravity. At early times, the enhancement of clustering 
is restrained since the Galileon self-interactions are the dominant ones suppressing 
their energy density relative to that of matter or radiation. Once the matter density 
has dropped sufficiently, the Galileon become an important contribution to the dark 
sector of the Universe driving cosmic expansion. 

The paper is organized as follows: We start in section 2 with a summary of the 
non-linear theory of massive gravity [27] which we use for constructing our proxy 
theory, where we discuss the successful implementation of ghost-free massive gravity 
in the decoupling limit. We then move to our proxy theory in section 3 where we 
first give some motivations on why we use a proxy theory rather than the exact 
non-linear theory. We then work out the consequences of this covariantization for 
cosmology and specifically for late-time acceleration in section 4 where we study 
also the stability conditions and the reason why a degravitating solution can not be 
constructed in our covariantized theory. We then move onto more general cosmology 
in section 5 and follow the helicity-0 mode contribution to the Universe throughout 
its evolution, before quickly presenting the consequences for structure formation in 
section 6. 
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2 Massive Gravity and its decoupling limit 



We use the same notation as [10], with (dn) 2 = g^d^nduTi, U^ u = D^Dyix and 
njL = g^n^II^ where the covariant derivative is taken w.r.t. g^ u and square 
brackets [...] represent the trace of a tensor [II 2 ] = IT^IP 1 ' and [II] 2 = n^IT^. 



2.1 Massive Gravity 

The first theory of massive gravity was proposed by Fierz and Pauli in 1939 [40], but 
was shown to be unstable by Boulware-Deser (BD) [41], due to the non-propagation 
of the Hamiltonian constraint at the non-linear level. This instability can also be 
seen in the Stiickelberg language [42] where in the decoupling limit, the helicity-0 
mode typically has equations of motion with more than two derivatives, and hence 
does not possess a well defined Cauchy surface. However it was shown in [27, 1] 
that the graviton potential could be built in such way as to remove any higher 
derivative term in the equations of motion, and obtain a Galileon-type of action for 
the additional helicity-0 mode. To review this, let us start with a graviton of mass 
m described by 

Mjj , ( „ m 2 



C = -fV=j\R- T U(g,H)j (1) 

with the potential U, where the tensor H^ u is constructed in terms of the metric 
g^ v and the four Stiickelberg fields <3> a by H^ v = g^ v — r] ab d IM ^ a d u ^ b . We can then 
split the Stiickelberg fields into helicity-1 and -0 contributions, but the helicity-1 
mode decouples in the decoupling limit and can be consistently set to zero. We 
therefore focus on the helicity-0 mode 7r and write $ a = (x a — ?7 aM <9 M 7r) such that 
= V + 2I V - V aP ^Jlp v . Defining the quantity K^g, H) = 5% — Vtf - Hi 
the most generic potential that bears no ghosts in the decoupling limit is 

U(g, H) = -4 (U 2 + a 3 U 3 + a 4 U 4 ) (2) 

where 0^4 are two free parameters and 

U 2 = [/C] 2 -[/C 2 ] (3) 
U 3 = [/C] 3 -3[/C][/C 2 ] + 2[/C 3 ] (4) 
U 4 = [/C] 4 -6[/C 2 ][/C] 2 + 8[/C 3 ][/C] +3[/C 2 ] 2 -6[/C 4 ]. (5) 

Notice that U4 can be expressed in terms of ^2,3 and the tadpole U\ = [K], [43]. 

It has then been shown that ghost-like pathologies in this theory of massive grav- 
ity theory disappear to all orders in the decoupling limit and at least up to quartic 
order beyond the decoupling limit, as well as completely non-linearly in some spe- 
cific cases. While this theory is completely covariant, studying the cosmology as 
well as other non-trivial curved geometries can be extremely complicated. Instead, 



5 



we will here focus on its version in the decoupling limit and covariantize the theory 
directly from this limit. We emphasize that the resulting proxy theory will be dis- 
tinct from the theory of massive gravity presented above, but presents nevertheless 
some interesting features for cosmology 

In the decoupling limit, taking the scales Mp\ — > oo and m — > 0, while keeping 
the strong coupling scale A 3 = M P im 2 fixed, the above Lagrangian (1) reduces to 
the more compact expression 

£ = —h^hap + h!»X$ + ^hrx$ + + i^f v T,u (6) 

where h^ u stands for the helicity-2 mode canonically normalized, is the Lich- 
nerowicz operator, the coefficients 02,3 are related to the free parameters 0^4 and 
xjn! 2 ' 3 ^ denote the interactions with the helicity-0 mode [27] 

xjS = DTT^-n^ (7) 

x$ = n^-nvrn^-l([n 2 ]-[n] 2 )^ (8) 
xf} = 6n^-6[n]n^ + 3([n] 2 -[n 2 ])n^-^([n] 3 -3[n 2 ][n] + 2[n 3 ]).(9) 

It is worth to mention that these interaction terms are all transverse and at most 
second order in time derivatives to ensure the absence of ghost. Being ghost-free, 
these interactions are closely related to the Galileon interactions and fulfill the same 
internal symmetry. In the next section we will covariantize these interaction terms 
h^X^J 2 and discuss their physical properties. 

In the case of external sources there is a coupling between the metric h^ v and the 
stress energy tensor via h^T^ but there is not such a direct coupling with the field 
7T. Nevertheless, if one diagonalizes the first interaction term h^xj^J by a change 
of variables of the form h^ u = h^ u + tttj^, then the coupling between ir and external 
sources will become transparent, ttT. 

3 Massive Gravity: A proxy theory 

Instead of studying the cosmology in the exact non- linear covariantized theory (1) 
which can be extremely hard, we use the alternative approach of covariantizing the 
Lagrangian in the decoupling limit, and use the resulting theory as a proxy. 

3.1 Covariantization 

We claim that the covariantized version of the above Lagrangian (6) is simply given 
by 

S = [ (Mpj_R + /7(vr, g,„) + £ matter (^, g^)) , (10) 
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where £ matter is the Lagrangian for the matter fields ip living on the geometry, the 
Lagrangian for 7r is 



C* = M P1 (-7TR - ^d^d^CT - ^d^TtU^L^ . (11) 
and the tensor L^ au ^ stands for the dual Riemann tensor 

jjwtvfi = 2 R^P + 2{R^g ua + R ua g^ - R» v g afi - R a/3 g^) + R{g^g aP - g^g va ) .(12) 

This form of tensor structure has been first discussed by Horndeski [44] in the context 
of the most general scalar-tensor theory, and more recently in [45, 46, 47]. However 
we point out here that these interactions come as a direct outcome of massive gravity. 
We can show explicitly the following correspondences, 

rig < — > -7iR (13) 
hTX® > -d^d^GT (14) 
hTX® <— > -d,nd u 7rU af} L^ . (15) 



The Einstein equation is given by 



G, v = M Pl T* u + T- atter (16) 

with 

and the structure of the Einstein and Riemann dual tensor ensure that 7r enters at 
most with two derivatives in the stress-energy tensor, 

rpir(l) -y-(l) , Q 

t;, (3) = x^ + h^n^^f (is) 

where are defined in (7 - 9). Furthermore, the fact that G 00 , G 0t , L° l0 ^ and 
j^oikj ^qtj-q a ^ most one time-derivative guarantees the propagation of constraints. 

Since we are not in the Einstein frame, these stress-energy tensor are only trans- 
verse on-shell, and satisfy the relation, D^T^ = d u n£ n where £ n is the equation of 
motion with respect to n. Since both the Einstein tensor and the Riemann dual 
tensor are transverse, this equation of motion is also at most second order in deriva- 
tive, 



£ 



5C n 

5tt 



-R - ^G^IV - ^-L^(IL^n a p + R\ ai/ d,nd^) = , (19) 
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where we have used the fact that 

D u D a D^L^ = -W pav d^L^ = i^TTAss (20) 

with £gb = R 2 + R 2 aai>p ~ ^Rfiv m the res ^ °f this paper, we study the result- 
ing cosmology in this proxy theory, starting with the existence of self-accelerating 
solutions. 



4 de Sitter solutions 

In what follows, we focus on the cosmology of the covariantized theory (10, 11), and 
focus for that on a FRW background with scale factor a(t) and Hubble parameter 
H . The resulting effective energy density and pressure for the field tt are then 

= M vi {QHt\ + QH 2 ti-^H 2 t\ 2 -^H^) (21) 



P n = 3Mi 



pi 



^H^{H + H 2 ) + ~#7f) + t^HHZH 2 + 2H) + AHn) 



-(7i(3H 2 + 2H) + 2Hn + n) , (22) 

and the equation of motion for tt (19) in the FRW space-time is equivalent to 

^ (3iJ 3 7r + 2HH-K + H 2 t) + ^ (3H 2 H tt 2 + + 2H 3 <kit) = R. (23) 

Similarly as in [13], this expression can be rewritten more compactly 

4> + 3H<p -R = (24) 

if we define <fi as 

*=*"(£* + ^"). (25) 
4.1 Self- accelerating solution 

Now, we would like to study the self- acceleration solution with H =const and H = 0. 
For the tt field we make the ansatz tc = q^-. Furthermore, we assume that we are 
in a regime where Hit <C tc so that we can neglect terms proportional to tt and 
consider only the terms including tt or tt. Thus, the Friedmann and field equations 
can be recast in 

2 

H 2 = — {6q - 9a 2 q 2 - 30a 3 <? 3 ) (26) 
3 

H 2 (18a 2 q + 5Aa 3 q 2 - 12) = 0. (27) 
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Assuming H ^ 0, the field equation then imposes, 



q = f 2 (28) 

6a 3 

while the Friedmann equation (26) sets the Hubble constant of the self-accelerated 
solution. Similar to what has been found in [2] our proxy theory admits a self- 
accelerated solution, with the Hubble parameter set by the graviton mass. For the 
stability condition of this self-accelerating solution the first constraint we have is 
to demand H > 0. The other constrain comes from the stability condition for 
perturbations on the background which we discuss in the following subsection. 

4.2 Stability conditions 

In the last subsection we have shown explicitly that our proxy theory exhibits a 
self- accelerating solution with H 2 w m 2 . Now, we would like to study whether the 
perturbations on this background are stable and what the constraints are. For this 
purpose, consider perturbations on the background solution of the following form 

7T = 7T (t) + &ir(t,x,y,z) (29) 

The second order action for the perturbations is 

£=-||V7r3^7rG^-|^^ (30) 
which can be written in the form 

C = K tt (5iT 2 -^(V57T) 2 ) (31) 
gt 



where 



3M Pl a 3 H 2 ( 6a 3 H . 
K * = Js K + ^-^J ( 32 ) 



01 = 1(2+ =r ). (33, 



and 

r 2 = x _ 

3 \ a 2 A 3 + Qa^H-k 

The condition for the stability is then given by K tt > 0, c 2 > and H 2 > 0, which 
are fulfilled if 

X 

a 2 > and > a 3 > --a\. (34) 

8 

To compare this result with the condition obtained in the decoupling limit, [2], we 
first mention that a\ eTe = — 2a2 here and a 3 ere = a3 here . In terms of the parameters 
used there, we need to compare our conditions a\ ete < 0, ag ere > —\{a^ vc ) 2 to the 
conditions af em < and -|(a| here ) 2 < af ere < -\{af ere ) 2 . We see that our theory 
is less constraining but is still within the parameter space derived in [2]. It is not 
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surprising that the stability condition in the decoupling limit and in our covariantized 
theory do not coincide totally as we have explicitly broken the symmetry when 
getting the proxy and in particular, our solution does now spontaneously break 
Lorentz invariance, which was not the case in the decoupling limit, [2]. 

It is also worth pointing out that the self-accelerating solution by itself does not 
propagate any superluminal mode, since 2/3 < (? s < 1. 

We emphasize as well that the constant tt solution is a dynamical attractor. For 
this we just consider time dependent perturbations 7r(i) = vr (t) + 8n{t) which is a 
special case of (29) fulfilling the same stability conditions. The equation of motion 
for perturbations simplifies to 

d t (a 3 5n) = (35) 

The solution for Sir is given by 

5ii(t)^a~ 3 . (36) 

Thus, these perturbations 5ir(t) redshift away exponentially compared to the tt = 
const self-accelerating solution. Therefore, the self-accelerating solution is an at- 
tractor. 

4.3 Degravitation 

More interestingly, one can wonder whether degravitation can be exhibited in these 
class of solutions. If one take tt = ir(t) and H = 0, it is straightforward to see that we 
obtain p n = 0, so the field has absolutely no effect and cannot help the background 
to degravitate. The situation is however different when setting tt = xo x a% a , which 
was possible in the decoupling limit. In this proxy theory, such a behaviour will 
involve explicit space dependances in the equations of motion which should be dealt 
with specifically. 

Interestingly, the interactions considered here are precisely of the same form as 
that studied recently in [46]. There as well, in the absence of spatial curvature k = 0, 
the contribution from the scalar field vanishes if H = 0. Comparing with [46], we 
can hence wonder whether the addition of spatial curvature k ^ in our proxy 
theory could help achieving degravitation, but relying strongly on spatial curvature 
brings concerns over instabilities which are beyond the scope of this study. 

5 Cosmology 

In the following we would like to discuss in more detail the interplay of all the 
constituents of the universe. We assume that matter, radiation and the scalar field 
7r contribute to the total energy density of the universe. 

H 2 = + p rad + p mat ) (37) 

3 
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Consider the scalar field ir as a perfect fluid with the effective energy density and 
pressure given by (21, 22). Thus, the equation of state parameter of this new field 
would be 

-l2a 3 H 3 n 3 + 2A 3 ij(2A 3 vr - a 2 n 2 ) + 2A 6 n + AHn(A 6 - 3a 3 Hn 2 - a 2 A 3 jr) 
Un ~ 3H(2AH - 10a 3 HH 3 + H(2AH - 3a 2 AH 2 )) 

3H 2 (2 A 6 tt - 7T 2 (a 2 A 3 + 6a 3 vf ) ) 
+ 3H(2AH - 10a 3 H 2 % 3 + H(2AH - 3a 2 AH 2 )) ' (38) 

At this point one should mention that the energy density for the 7r-field is not 
conserved but rather given by D^Tjf = d u 7t£ n (where £ n is the equation of motion 
for 7r), which is not surprising since 7r is non-minimally coupled to gravity in the 
Jordan frame. Therefore, we can have w 7r < — 1. 

In the following we will first assume that at early times in the evolution history of 
the Universe we can neglect the extra density coming from the helicity-0 p n . We 
will then check this assumption by plugging the solution for H back in the equation 
of motion for n. If we assume that at early times the radiation density dominates, 
we simply have 

H 2 = ^G_ p r,A a -A a ^ t l/2 w = l/ g (3g) 

3 

During the radiation era, the dominant terms in the equation of motion for 7r are 
then ^H 2 Hn 2 + ^H 4 tt 2 + ^HHn = which can be solved assuming the 
previous expression for H (39) 

TW-t 1 ' 75 yielding /O-MpdT 1 / 4 (40) 



At later times when the matter dominated epoch starts we have 

#2 = ^ p mat a -3 fl „ f 2/3 CO = (41) 

3 

Now the dominant terms in the equation of motion for 7r are ^pp-if 3 7T + -j^HHjt + 
^-H 2 n - 12H 2 - QH = 0. We get for vr this time 

i 2 A 3 w i 3M Pl (-Ua 2 + 5a 3 )A 3 lA , 

7r mat ~ c 2 ■ t + — yielding = c^t' 1 + m 2 ^— . 42 

4a 2 32a 2 

Summarizing, during radiation domination the effective energy density for the 7r-field 
goes like p^ ad ~ t -1 / 4 while during matter domination as p^ at ~ t^ 1 and approaches 
a constant at late time. As shown in the figure p 71 can be neglected at early times 
where p^ at < p mat and p^ d < p rad . 
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Figure 1: Fluid densities p rad ~ a -4 , p mat ~ a -3 and p n during the epochs of ra- 
diation, matter and A-domination normalised to today p^. During the radiation 
domination the energy density for ir goes as p^ ad ~ a" 1 / 2 and during matter domi- 
nations as p^ at ~ a~ z l 2 and is constant for later times p\ = const. 



6 Structure formation 

We end the cosmological analysis by looking at the evolution of matter density 
perturbations. The density perturbations follow the evolution 

5 m + 2H5 m = — f (43) 

where ip is the Newtonian potential. The effects of 7r are all encoded in its contribu- 
tion to the Poisson equation. Consider perturbations of ir around its cosmological 
background solution n(x,t) = n (t) + <f>(x,t). gives a contribution to the Newto- 
nian potential of the form if) = 0/Mpj. In the Newtonian approximation we have 
|0| <C |V0|. The equation of motion for the scalar field in first order in is 

- ^GTD^ - ^-L^AUfD^D^ + 2R^ d^d^) = SR (44) 

which is equivalent to (neglecting 0) 

^ = SR (45) 



^(3H 2 + 2H) + ^(2HH + 2HHtt + HH ) 
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and last but not least we need the trace of Einstein equation, (16-18). Perturbing 
the trace to first order, we get 



M Pl 6R 



a 2 



Q 

3 - 2-^|(2#7r + 7f ) - ^(2^7r 2 + 5HV + 4Hn n ) 
5T 

+ Afe (46) 



To reach that point, we have neglected the perturbations of the curvature of the 
form SRttq as they are negligible compared to Mp\5R since we work in the regime 
where ttq <C Mp\- We have also ignored terms of the form 5 R(dito) 2 / A 3 , which 
could a priori be relevant, but their inclusion would require solving the full Einstein 
equation, which is beyond the scope of this study. As a first approximation, such 
terms are hence ignored. 

The perturbations for the source is just given by 5T = —p m $m for non-relativistic 
sources, thus we have 

V 2 p m S m 



a 2 ~ 3M Pl Q 

where Q stands for 



(47) 



Q = l--r£(2Hir + it + M Pl (2H + 3IP)) (4* 
A d 

°3 ( r r T 2-2 , o ir--2 16M P1 M r/3 • , o U H ■ , zj2 ■ 



v . . 5#X + 2H Pio + 4F7f 7r ^(2H s tt + 2HHir + H 2 jf Q ] 

Finally, the modified evolution equation for density perturbations is 

L + 2H5 m = ^(l + ^] ■ (19) 



Knowing the background configuration it is then relatively straightforward to derive 
the effect on structure formation. We recover the usual result that when the field is 
screened, Htt > A 3 , the extra force coming from the helicity-0 is negligible and the 
formation of structure is similar as in ACDM. 



7 Covariantization from the Einstein Frame 

After having studied the cosmology and the structure formation of our proxy theory 
in the Jordan frame, the natural question is whether we would expect similar results 
in a different frame. Instead of covariantizing our Lagrangian in the Jordan frame, 
it is on an equal footing to go to the Einstein frame where the Ricci scalar is not 
multiplied by the scalar field 7r and covariantize the theory at that stage. Since it is 
unclear at first sight which frame is the physical frame we will consider both frames 
and discuss their differences. Our starting Lagrangian was 

c = -^effhap + rxg + + ^rxg + ^-^^ (so) 
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Now, when we do the following change of variables 

hfw = h^u + nr]^ (51) 

we can diagonalize the first mixed term h^Xfi) such that the Lagrangian takes the 
following form 



+^X<£ - 2^(<9vr) 2 ([n 2 ] - Dvr 2 ) + _( V + (52) 

Covariantizing this action is straightforward. We use again the correspondences in 
(3.1) and it has been shown explicitly that the covariant equivalence to nir(dn) 2 
and -2(<97r) 2 ([n 2 ] - Dtt 2 ) are given by D7i(dn) 2 and 2(dn) 2 (Dn 2 - [II 2 ] - \{d-K) 2 R) 
respectively, which do not yield any ghostlike instabilities ([23], [11]). Thus, the 
covariantized action in the Einstein frame is simply given by 



C = M 2 Pl R+ 3 -nnn-^d,nd^- 3 -^nn(dn) 2 

-^d^Il^L^ + 2^(dn) 2 (Dn 2 - [II 2 ] - \{d*fR) (53) 
+£ m (V>, (1 + vr)fiv) 

Similarly as before, the properties of G^ v and L^ au ^ ensure that their equations of 
motion lead at most to second order derivative terms. To find a self-accelerating 
solution we set again a pure de Sitter metric, with tt = qA 3 /H. The Friedmann and 
the field equations then take the form 

Mp\H 2 = 3M Pl a 2 A 3 q 2 + 10M P1 a 3 A 3 g 3 - + 3a 2 -^g 3 + 15a 3 -^g 4 

A 3 (-l + 3q(a 2 + 4a 3 g)) + 2M P1 (a 2 + 3a 3 q)H 2 = (54) 

When comparing the above Friedmann and the field equations with the one we had 
in the Jordan frame (26), we see significant differences coming from the extra terms 
which were not there in the Jordan frame. These terms yield Friedmann and field 
equations proportional to g 4 and if 4 which are more difficult to solve. 
For fairness, we should compare both actions in the same frame. We do so by 
performing a conformal transformation on the action (11): 

= Vfg^ with Q 2 = ( 1 - J . (55) 



pi 



For simplicity we consider the case for which a 3 = 0, so under this conformal trans- 
formation the covariantized action (11) becomes 

Cj = M 2 ,R - 3 -Q-\dn) 2 - ^ (WttG,, + —{^Dk + 5^(^)^(56) 
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In the limit where ir <C Mp\ we have then finally the following expression 

Cj = M&R + 3 -nDn - ^d,nd u nG^ - \%^{d-nf , (57) 

which coincides with the theory obtained from the Einstein frame, (53) with 03 = 0. 
So within the regime of validity of our results, our conclusions are independent of the 
choice of frame. However beyond this regime the theory originally constructed from 
the Jordan frame could violate the null energy condition from the term proportional 
to (<9tt) 4 . 

8 Discussion 

Common to generic massive gravity models is the presence of at least one addi- 
tional helicity-0 degree of freedom which originates from a spin-two graviton field. 
In contrast to an arbitrary scalar field, its dynamics is not driven by a potential but 
rather by specific derivative terms fixed by symmetries. Such a degree of freedom 
is experimentally testable. On solar system and galactic scales gravity is very well 
compatible with GR and correctional terms from a modified gravity model are ex- 
cluded. On these scales, however, the effects of massive gravity might be cloaked 
by the Vainshtein mechanism, where the helicity-0 mode interactions become appre- 
ciable to freeze out the field fluctuations. A number of screening mechanisms have 
been devised such as the chameleon and symmetron mechanisms. Contrary to those 
the success of the Vainshtein mechanism does not rely on some specific potential 
but instead on derivative interactions, which cause the helicity-0 mode to become 
decoupled from matter and light on short distances, which nevertheless could have 
observational signatures on larger scales in cosmic structure formation. 

A further benefit of the Vainshtein mechanism would be a natural solution to 
the problem that the matter density and the cosmological constant are of similar 
magnitudes today. The cosmological application of the Vainshtein mechanism works 
such that at early times, the galileon interactions are dominated by self-interaction, 
which suppresses their energy density relative to that of matter or radiation. If the 
matter density has decreased sufficiently by cosmic expansion, the galileon consti- 
tutes an important contribution to the energy density of the universe and drives 
cosmic expansion. 

Although the precise cosmology of the theory of massive gravity proposed in 
[1] has not been derived, most models of massive gravity give rise to alterations of 
the Hubble function similar to DGP-gravity. Experimental tests of the expansion 
dynamics of the Universe include the distance-redshift relation of supernovae, and 
measurements of the angular diameter distance as a function of redshift, as in the 
case of the cosmic microwave background and the baryon acoustic oscillations. Apart 
from these geometrical tests, the time-dependence of evolving cosmic structures 
can be investigated, and the influence of the gravitational theory on the geodesies 
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of relativistic (photons) and nonrelativistic (dark matter) test particles. The first 
category includes gravitational lensing and the Sachs- Wolfe effects, which have been 
shown to differ from their GR-expectation in some modified gravity theories, and 
a similar result can be expected in massive gravity. The second category includes 
the homogeneous growth of the cosmic structure, and the formation of galaxies and 
clusters of galaxies by gravitational collapse. Again, massive gravity would influence 
the time sequence of gravitational clustering and the evolution of peculiar velocities, 
as well as the number density of collapsed objects. In particular, massive gravity 
would enhance gravitational clustering since they tend to lower the collapse threshold 
for density fluctuations in the large-scale structure, leading to a higher comoving 
number density of galaxies and clusters of galaxies. Naturally, these changes are 
degenerate with a different choice of cosmological parameters and with introducing 
non-Gaussian initial conditions, which would be very interesting to quantify. Recent 
discrepancies of ACDM with observational data on large scales include the number of 
very massive clusters, the strong lensing cross section, anomalous multipole moments 
of the CMB, the axis of evil, and large peculiar velocities. It is beyond the scope 
of this paper to address these issues using massive gravity, but we propose how 
to proceed from constructing a massive gravity theory to providing observationally 
testable quantities. From our point of view it is advisable to focus on probes of large 
scales, due to difficulties related to nonlinear structure formation and the influence 
of baryons on small scales. Natural questions concern the homogeneous dynamics of 
the Universe, the formation of structures and the shape of geodesies of relativistic 
and nonrelativistic particles. Basically, using our proxy theory one should be able 
to make predictions concerning these four issues and the combination of the four 
should give insight into the nature of the gravitational sector. In our proxy theory, 
from the modified field equation the Hubble-function can be derived easily, which 
allows the definition of distance measures, needed in the interpretation of supernova 
data. Cosmic structure growth tests the Newtonian limit for slowly-moving particles 
and describes the clustering of galaxies and the growth of structures which are 
investigated by e.g. gravitational light deflection. Lensing, in turn, makes use of 
the geodesic equation for relativistic particles, and measures the correlation function 
of the matter density, weighted with the lensing efficiency function, which in turn 
is derived from distance measures. It would be quite interesting to study these 
observational consequences and constraints of our proxy theory in a future work. 
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